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ABSTRACT 

Chandra observations of NGC 4636 show disturbances in the galaxy X-ray halo, including arm-like high 
surface brightness features (tentatively identified as AGN driven shocks) and a possible cavity on the west side 
of the galaxy core. We present Chandra and XMM-Newton spectral maps of NGC 4636 which confirm the 
presence of the cavity and show it to be bounded by the arm features. The maps also reveal a ~15 kpc wide 
plume of low temperature, high abundance gas extending 25-30 kpc to the southwest of the galaxy. The cavity 
appears to be embedded in this plume, and we interpret the structure as being entrained gas drawn out of the 
galaxy core during previous episodes of AGN activity. The end of the plume is marked by a well defined edge, 
with significant falls in surface brightness, temperature and abundance, indicating a boundary between galaxy 
and group/cluster gas. This may be evidence that as well as preventing gas cooling through direct heating, AGN 
outbursts can produce significant gas mixing, disturbing the temperature structure of the halo and transporting 
metals out from the galaxy into the surrounding intra-group medium. 

Subject headings: galaxies; individual: NGC 4636 — galaxies; intergalactic medium — galaxies; halos — 
X-rays: galaxies 


1. INTRODUCTION 

NGC 4636, the dominant galaxy of a group on the out¬ 
skirts of the Virgo cluster, is one of the most well-known 
nearby ellipticals. It is also one of the most X-ray lumi¬ 
nous, and has been studied in detail by every major X- 
ray imaging observatory to date. Spectral imaging with 
Einstein showed the galaxy to be surrounded by an ex¬ 
tensive halo of hot gas with a temperatur e of 0.78-1.21 
keV (90% eiTor bounds, lEorman et al.lll985]) . Eurther stud¬ 
ies using ROSAT and ASCA measured the metal abun¬ 
dance of the gas, found gradients in both the tempera¬ 
ture and abundan ce with radius, and showed the halo to 
be very extended (|A^aMeLalJ[19^ |Tlinchi^i_eL^ll994t 
Matsus iita et alJ 119971 iFmoguenov&JonesI I2000t iBuotel 
2000aHO’SullivanetalJl^03al) . More recently. NGC 4636 
has been observed by both Chandra and XMM-Newton. The 
initial Chandra ACIS-S observation showed unusual struc¬ 
tures in the core of the galaxy, most notably “spiral arm” 
features and a possible ring of high surface brightness emis¬ 
sion, thought to be the product of shocks driven by a previous 
AGN outburst (llones et al.ll2002ft . Detailed spectral analysis 
of these data suggests the presenc e of a cavity in t he X-ray 
halo on the west side of the core JOhto et al.ll2003l) . Again, 
AGN activity is the likely cause; plasma from the radio jets 
of the AGN can displace the gas of the X-ray halo, leav¬ 
ing an apparent void and altering the projected X-ray prop¬ 
erties at that point. Analysis of XMM RGS spectra shows 
that there is little or no emission from gas at temperatures 
belo w 0.5 keV, and therefore no cooling flow in the sys¬ 
tem (IXu et al.ll2002h . NGC 4636 hosts a weak radio source 
(1.4x10^* ergs”') with small scale j ets which account for 
~40% of the radio emission at 1 .4 Ghz jBirkinshaw & DaviesI 
Il985t[^nger & WarwiclJl98^. The Ha lum inosity is mea¬ 
sured to be 2.3 X10^^ erg s~' JHo et al.lll997h . and an upper 
limit of 2.7x lO*'* erg s“* is estimated for the X-ray emission 
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of the AGN (iLoewenstein et alJ2001h . 

Using the high quality data available from both Chandra 
and XMM, we present spectral maps of the inner galaxy halo 
of NGC 4636, which reveal interesting temperature and abun¬ 
dance structures. We consider the origin of these structures 
and their impact on our understanding of galaxy and group 
halos. 

2. OBSERVATIONS AND DATA ANALYSIS 

NGC 4636 was observed with the Chandra ACIS-S (chip 
S3) for 53 ks on 2000 January 26-27 (ObsID 323) and 
with XMM EPIC for 64 ks on 2001 January 5-6 (ObsID 
0111190701). Eor the Chandra observation, the instrument 
operated in faint mode, and reduction and analysis of the 
data were performed using CIAO v.3.1 and CALDB v.2.28. 
Reprocessing, removal of background flares, and creation of 
responses and background w ere carried out as described in 
lO’ Sullivan & PonmanI ll2004ll . Point sources were identified 
using the WAVDETECT tool. Apart from those in the galaxy 
core, which were considered false, all were removed. 

The XMM EPIC instruments were operated in full frame 
mode with the medium optical Alter. The raw data were 
processed using SAS v.6.0, and reprocessing. Altering for 
bad pixels and columns and times of high background, 
and point source remov al were performed as described in 
lO’ Sullivan et al.l j2003bl) . including use of the ‘(FLAG == 
0)’ Alter. Background events lists were generated from the 
iRead & PonmanI (I2003h blank-sky background data, cleaned 
to match the data. The SAS EVIGWEIGHT task was applied 
to both source and background events lists, “correcting” the 
events to account for vignetting. This allows the use of sin¬ 
gle, on-axis ARP response files for each instrument, which 
were generated using ARFGEN. Pregenerated RMP files were 
used, selected by position. The exposure times after cleaning 
were ^58.4 ks for the MOS cameras, ~50.9 ks for the PN, 
and ~38.9ks for ACIS-S. 

In order to look for correlations between the structures vis¬ 
ible in the X-ray surface brightness images and temperature 
and abundance changes in the halo, we created spectral maps 
from the Chandra and XMM datasets using the method de- 
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Fig. 1.— Upper panel: XMM spectral maps of the southwest part of the halo of NGC 4636, showing (left) Temperature in keV and (right) Abundance in 
solar units. The maps use 75x73 8.8'^^ pixels, whose spectral extaction regions are squares of side 9-67.5^^^. The diamonds on each map mark the center of the 
galaxy, regions marked 1 and 2 are, respectively, the cavity and plume referred to in the text. Black lines mark the region of overlap with the Chandra spectral 
maps (see Figure|3- Lower panel: En'or maps showing the 90% upper and lower limits on temperature and abundance in each pixel, and a gaussian smoothed 
MOS 1 surface brightness image of the same field. Lines mark the southern surface brightness edge referred to in the text. 10 kpc is equivalent to ~130". 


scribed in lO’Sullivan et all j2005l) . The spectra for each map 
pixel were required to have >800 counts, and were fitted with 
an absorbed MEKAL model with hydrogen column fixed at 
the galactic value (1.8 x 10^** cm“^). Energies below 0.35 (0.5) 
keV and above 6.0 (8.0) keV were ignored in the Chandra 
{XMM) fits. Pixels whose 90% errors on kT were >10% 
(15%) were excluded from the Chandra {XMM) maps. Eig- 
ure^ shows temperature and metal abundance maps gener¬ 
ated from the XMM data, with associated 90% error maps. 
Eor comparison, Eigure |2l shows Chandra temperature and 
abundance maps, plotted using the same color scale as used 
in EigureH] as well as an X-ray surface brightness image. In 
each hgure, all panels are plotted on the same scale and coor¬ 
dinates, so as to be directly comparable. 

Several important features are visible in the temperature 
maps. Perhaps the most obvious is a region of low tempera¬ 
tures (blue-purple) in the galaxy core and extending in a curve 
to the south and southwest. To the west of the galaxy core is 
an area of higher temperatures apparently surrounded by cool 
gas (region 1 in Eigure^. To the west of this is a large area 
of moderately cool temperatures (region 2) apparently con¬ 
nected to the core on both sides of region 1. To the north, east, 
and southeast of the core the temperatures are higher, and this 
high temperature gas appears to form a boundary around the 
cool galaxy core and regions 1 and 2, though at the western 
edge of region 2 it is appears to be relatively narrow, with 
cooler gas beyond. Three main features can be drawn from 
the abundance maps; first that there is a boundary between 
high abundance gas in the inner halo, and low abundance gas 
at the southern and western edges of the maps, with a similar 
boundary between high and low temperatures; second that the 


highest abundances are seen to the southwest of the galaxy 
core, near region 2 of the temperature map; third, that the 
abundances in the galaxy core appear to be only ^^0.5 solar. 

Comparing these features to the X-ray surface brightness 
images we hnd that the high temperature region 1 seems to be 
bounded by the western part of the “spiral arm” features. In 
particular, the southwest “spiral arm” marks the boundary be¬ 
tween region 1 and the low temperature curve extending from 
the galaxy core (see Eigure |2j. On larger scales the surface 
brightness appears roughly circular (see Eigurelfl, but a sharp 
drop in surface brightness along the south and southwest sides 
of the inner halo is visible, coincident with the end of region 
2 and the southern edge of the area of high abundance. This 
is probably the boundary between the galaxy halo and a sur¬ 
rounding intra-group or -cluster medium (ICM). Regions 1 
and 2 lie roughly along the minor axis of the stellar body of 
the galaxy. The AGN radio jets do not extend as far as region 
1, and are i nitially aligned ^43° away fr om the minor axis of 
the galaxy llBirkinshaw & Dayieslll98-^ . but at a small dis¬ 
tance from the center they appear to bend, with the southern 
jet turning to point toward region 1. 

There are issues which could affect the absolute accuracy 
of our maps. Many map pixels will not be independent, due 
to oyerlapping spectral extraction regions and the large XMM 
PSE. We haye also not corrected for differences in the soft 
component of the background between the source and back¬ 
ground datasets. The fact that both Chandra and XMM maps 
show the same features and similar temperatures and abun¬ 
dances at giyen positions is a strong argument that the PSE 
issue and background differences are not causing signihcant 
problems. It is also important to point out that we are in- 
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Fig. 2.— Chandra temperature (A) and abundance (B) maps of the core of NGC 4636, using the same color scheme as in Fig.0 The maps use 64x64 6.4"^ 
pixels, whose spectral extaction regions are squares of side 8.4-63”^. For comparison, a Chandra 0.5-2.0 keV surface brightness (C) image is also shown. All 
images are on the same scale and centered on the same coordinates. VLA 21cm radio emission contours are overlaid on panel C, and outlines of the surface 
brightness features are overlaid on panels A and B to aid comparison. 10 kpc is equivalent to '--'130''. 


terested in relative differences in temperature and abundance 
across the maps, rather than the absolute value in each pixel. 
However, as a further test, we carried out more traditional 
spectral hts for selected regions of the XMM data. We ex¬ 
tracted spectra from regions centered on the features in the 
spectral maps and fitted them with a variety of models. With 
the exception of the galaxy core, all hts produced results in 
agreement with the maps. The galaxy core required a multi¬ 
temperature plasma and power-law model, and while the re¬ 
sulting temperatures agreed with the maps, abundance was a 
factor ~2 higher. The cause is probably the “Fe-bias” effect, 
the single temperature map ht u nderestimating the abundance 
of mul ti-temperature emission (iBuote & Fabianlfl998t iBuotel 
l2000bl) . With this exception, we believe these spectral hts 
show the maps to be an accurate representation of the true 
temperature and abundance structure. 

3. DISCUSSION 

The high temperatures observed at region 1 corre spond 
closely with the results reported by lOhto et aTI ll2003l) . and 
conhrm their suggestion of a cavity in the gas. The apparent 
curve of the southern radio jet toward region 1 adds further 
support, and we note that the cavity lies along the minor axis 
of the galaxy, where the density gradient might be assumed 
to be steepest. The temperature map suggests that the cav¬ 
ity is non-spherical, and instead has a concave surface toward 
the galaxy center. This is reminiscent of the bell-like shapes 
seen in some numerical simulations of buoya ntly rising bub¬ 
bles o f radio plasma in galaxy clusters (e.g.- IChiirazov et^ 
120011) . Giyen the high temperatures around the north, east, 
and southeast of the core, it seems likely that the core is sur¬ 
rounded by a spherical shell of hot gas through which a plume 
of cool gas, ending in region 2, has penetrated. Under this as¬ 
sumption, the high temperature of the cayity is a projection 
effect, indicating that there is less cool gas along the line of 
sight than is the case for the surrounding regions. 

Taking a simple model of the cayity as a spherical bubble 
of radio plasma expanding into and rising through the X-ray 
plasma of the galaxy halo, we can estimate the time taken 
to create the cayity. Assuming the cayity expands without 
causing shocks in the halo, we can limit the expansion ye- 
locity to the sound speed of the halo gas, leading to an ex¬ 
pansion time of ~23 Myr. From l^iirazov et al.l (l20nnl) . the 
yelocity at which the bubble rises owing to its own buoyancy 
is Vb = C\JrGM/R^, where R is the distance of the bubble 
from the galaxy core, r is the radius of the cayity, and C 
is a numerical constant for which a yalue of 0.5 is accepted 


for incompressible fluids. This yields a buoyant rise time 
of '^47.5 Myr, but this should be considered a lower limit 
as this assumes a large difference in gas density between the 
galaxy core and the current location of the bubble. These two 
timescales together suggest that the bubble may haye formed 
close to its current location, rather than at the core of the 
galaxy. We can also estimate the power required to produce 
the cayity, assuming it expands at the sound speed and the 
pressure equilibri um between the radio a nd X-ray plasmas is 
maintained. From IChurazoy et°^ (l200(]l) . the power required 
is Power = KP/ht where X = (7- 1)/7, 7 is the adiabatic 
index of the relatiyistic gas in the bubble (i.e. 7=4/3), t is the 
expansion time of the cayity and P is the pressure of the halo 
around the cayity. From this we estimate the power required 
to be ~4.8xl0"^^ erg s“'. This can be compared to the cur¬ 
rent luminosity of the core of NGC 4636, Lx—3 x 10'*' erg s“' 
within 20 kpc. 

These estimates are considerably higher than the obseryed 
X-ray and radio emission of the AGN, but 2-3 orders of mag¬ 
nitude below estimates of the power requir ed to create the 
larger cayities seen in some galaxy clusters (iBohringer et alJ 
Hool. The radio luminosit y to jet power conyersion fac¬ 
tors of iBicknell et all Jl997h suggest that the total current 
jet power is < 3 x 10** erg s“' (assuming a magnetic held 
strength >10“*G and the age of the jet to be the bubble ex¬ 
pansion timescale ), and giyen the lack of a direct connection 
between the jets and the cayity, it seems likely that the AGN 
is currently in a quiescent state, and that the cayity is a rem¬ 
nant of a preyious outburst. A similar situation is obseryed 
in Abell 4059, which has large cayities in the X-ray halo, but 
a ce ntral radio galaxy w hose has declined since their forma¬ 
tion dHeinz et alJl2002ft . There is also the possibility that the 
AGN outburst was more powerful than our calculations sug¬ 
gest, and caused supersonic expa nsion of the bubble ; the “spi¬ 
ral ann^featoesjr^e core were LTones et alJll2002h identihed 
by lTones et alJ J2002ft as shocks, and a rapidly expanding cay¬ 
ity might shock surrounding gas. Howeyer, the resolution of 
our spectral maps and the size of the spectral extraction re¬ 
gions means that we cannot comment on these features, ex¬ 
cept to point out their association with the cayity. In any case, 
the cayity size and power required are small by the standards 
of galaxy clusters. 

One possible reason for this could be that AGN power 
is related to the mass of the central supermassiye black 
hole (SMBH), so that differences in SMBH mass (Msmbh) 
cause differences in ou tburst strength and cayity size. 
iMerritt & Ferrare^ J2001I) estimate Msmbh=7.9 x 10^ Mq for 
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NGC 4636, based o n the stellar velocity dispersion, while 
Ivan der Marell lll999ll calculates Msmbh=3.6x 10* Mq based 
on the central optical brightness p rofile. Estimates for M 87 
are higher, Msmbh=3.0x 10^ Mq JTremaine et alJlT002h be¬ 
ing typical. This suggests that NGC 4636 has a SMBH ~5- 
40 times less massive than M87, and if outburst power is di¬ 
rectly proportional to Msmbh, this might explain the differ¬ 
ence in cavity sizes. It is unclear whether there is a general 
difference in Msmbh between group and cluster dominant el¬ 
lipticals, as only a few measurements are available . How¬ 
ever, we note that in the iMerritt & Ferrare^ j2001h sample, 
all galaxies with Msmbh >10^ Mq are cluster dominant el¬ 
lipticals, apart from NGC 4649, which is at the center of a 
subclump of the Virgo cluster. Another point is that galaxy 
groups (and individual ellipticals) have short cooling times 
compared to clusters. This implies that the AGN duty cy¬ 
cle in NGC 4636 should be short; after an outburst a cooling 
flow will be rapidly reestablished and may soon fuel a new 
outburst. These factors suggest a model in which NGC 4636 
produces small outbursts relatively often, while cluster dom¬ 
inant ellipticals produce a smaller number of more powerful 
events. 

Region 2 may be relevant to this argument. The temper¬ 
ature maps suggest that there is a plume of gas extending 
southwest from the galaxy core, with region 2 forming the 
tip of the plume and the cavity of region 1 embedded within 
it. The direction of this plume matches that of the galaxy mi¬ 
nor axis. The high abundances to the southwest of the galaxy 
core suggest that much of the gas in this area has been en¬ 
riched, and the most likely explanation of this is that the gas 
was enriched inside the galaxy, then transported outward. In 
such a model, the plume and cavity can be interpreted as a 
snapshot of an ongoing series of small AGN outbursts, pro¬ 
ducing cavities which entrain metal rich, cool gas from the 
galaxy core. Entrainme nt of gas by buoyant b ubbles is pre¬ 
dicted from simulations dChiirazov et al.ll^0flh . and high res¬ 
olution X-ray observations of the radio/X-ray arms of M87 
show regions of multi-temperatu re gas which seem likely to 
be entrained gas llMolendill2002h . The width of these struc¬ 
tures in M87 is ^8-10 kpc, whereas the NGC 4636 plume is 
~15 kpc across. 

If this scenario is correct, the cool plume and high abun¬ 
dances require the AGN of NGC 4636 to have been active 
prior to the episode which produced the cavity visible in the 
temperature map, and to have significantly affected the struc¬ 
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ture of the galaxy halo. The removal of cool, enriched gas 
from the galaxy core would clearly affect the formation of 
any cooling flow, and could have important implications for 
enrichment in galaxy groups, as it provides a mechanism for 
the diffusion of metals from the central elliptical into the sur¬ 
rounding intra-group medium (IGM). The sharp drop in X-ray 
surface brightness, temperature and abundance at the south 
and west edges of the maps all suggest that this marks the 
boundary between gas associated with the galaxy, enriched 
by stellar mass loss and supernovae, and gas associated with 
a larger potential, which has a much lower metallicity. This 
surrounding gas could be the halo of the group of which NGC 
4636 is the dominant member, or the Virgo ICM. Confirma¬ 
tion of ICM with kT~0.5 keV and O.1-O.2Z0 abundance at 
2.6 Mpc from the cluster core would be an interesting re¬ 
sult in itself, but may be unrealistic. Simulations of entrain¬ 
ment by AGN jets in clusters suggest that the mixing effect 
is not strong en ough to affect observed abundance gradients 
( {Bijjggeni2flfl2h . However, these simulations do suggest that 
gas can be moved distances of 15-20 kpc, comparable to the 
motion suggested by our maps. If there is a boundary in 
NGC 4636 between gas associated with the galaxy and the 
IGM at a distance of 25-30 kpc, then motions on the scale 
seen in the simulations would be sufficient to mix the two, 
particularly if multiple AGN outbursts are considered. 

4. CONCLUSIONS 

Using X-ray data from XMM-Newton and Chandra, we 
have mapped the temperature and metallicity structure of the 
inner halo of NGC 4636. The maps confirm the presence 
of a cavity to the west of the galaxy, and show a plume of 
cool, metal-rich gas extending beyond the cavity to the south¬ 
west. Both cavity and plume appear to be the product of 
past AGN activity, the AGN being quiescent at present. The 
most likely scenario involves AGN outbursts producing bub¬ 
ble which then entrain enriched gas as they rise buoyantly 
from the galaxy core. If this is the case, then AGN activity 
may play an important role in producing both the temperature 
and abundance structures we observe in galaxy groups. 
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